Disentangling Scaling Properties in Anisotropic Fracture 
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Structure functions of rough fracture surfaces in isotropic materials exhibit complicated scaling 
properties due to the broken isotropy in the fracture plane generated by a preferred propagation 
direction. Decomposing the structure functions into the even order irreducible representations of 
the SO(2) symmetry group (indexed by m = 0,2,4...) results in a lucid and quickly convergent 
description. The scaling exponent of the isotropic sector (m = 0) dominates at small length scales. 
One can reconstruct the anisotropic structure functions using only the isotropic and the first non 
vanishing anisotropic sector (m = 2) (or at most the next one (m = 4)). The scaling exponent of 
the isotropic sector should be observed in a proposed, yet unperformed, experiment. 



Introduction: Imagine a fracture experiment in 
which an initial circular cavity is made to propagate by 
a tensile load such that the crack edge remains circular 
on the average, without any preferred propagation direc- 
tion in the plane normal to the load, see Fig. ^ From 
the point of view of the scaling properties of the rough 
fracture surface that is left behind the advancing crack, 
such an experiment is the analogue of homogenous and 
isotropic turbulence in nonlinear fluid mechanics. Nor- 
mal experiments in both fracture and turbulence involve 
symmetry breaking; the boundary conditions introduce 
anisotropy, making the discussion of scaling properties 
non-trivial. In turbulence it was shown how to disen- 
tangle the anisotropic contributions from the isotropic 
one by projecting the measured correlation and struc- 
ture functions on the irreducible representations of the 
SO (3) symmetry group IJ. The scaling phenomena seen 
in the isotropic sector of anisotropic experiments are 
identical to those expected in the hypothetical experi- 
ment of homogenous and isotropic turbulence. The aim 
of this Letter is to introduce a similar concept to the field 
of fracture: we will show that decomposing the height- 
height structure functions of fracture surfaces into the 
irreducible representations of the SO (2) symmetry group 
results in a simplification and rationalization of the scal- 
ing properties that is not totally dissimilar to the one 
obtained in turbulence. The scaling properties of the 
isotropic sector should be observable in principle in an 
experiment like the one shown in Fig. ^ which contrary 
to turbulence may be performed in reality. 

The pioneering experimental work described in Ref. 
PI drew attention to the fact that fracture surfaces are 
graphs in 2+1 dimensions when the broken sample is 
three dimensional. This initial insight was followed up 
by a considerable number of works ^^J that focused on 
the scaling properties of such graphs under affine trans- 
formations. In 2-f 1 dimensions one denotes the graph as 
h{r) and considers the structure function 5*2 (•^), 



S2{e)^mr+l)-h{T)f) 
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where angular brackets denote an average over all r. 
Initially no attention was paid to the fact that the 



FIG. 1: Sketch of a hypothetical fracture experiments ar- 
ranged to allow a crack to develop in an isotropic fashion, i.e. 
with all angles 6 being statistically equivalent. On the left, the 
full three dimensional experiment is shown, where the tensile 
axis is along z and a circular cavity is in the xy plane. On 
the right, a magnified version of the circular cavity in the xy 
plane is shown. 



isotropy in the fracture plane is broken due to initial 
conditions that lead to a preferred propagation direc- 
tion, and the statement of was that the structure 
function is a homogeneous function of its arguments, 
5*2 (A£) ^ A''*^'S'2(€) . In fact such a statement is ten- 
able only if the fracture process and the material are 
isotropic. Usually the crack propagates predominantly 
in one direction (say x) and the vector £ defines an angle 
9 with respect io x, 9 — cos~^(& • t). There is no reason 
why the scaling exponent C,'^'^\ if it exists at all, should 
not depend on the angle 9. Indeed, in the later work 
that followed ,2] this problem was recognized and scaling 
exponents were sought for one dimensional cuts through 
S2{i): typically parallel and orthogonal to the direction 
of the crack propagation. Besides the obvious meaning of 
'parallel' and 'orthogonal' to x, no reason was ever given 
why these particular directions are expected to provide 
clean scaling properties. We argue below that in general 
such one dimensional cuts exhibit a mixture of scaling 
exponents with amplitudes that depend on the angle ^, 
where 9 = Q and 9 — it/ 2 are not special. 
Approach: Given an experimental surface h(r) we first 
compute the second order structure function defined in 
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Eq. (Q. The vector i is associated with a norm ^ and 
an angle Q. By construction, the second order structure 
function is symmetric under d — > + tt. Accordingly, 
decomposing the structure functions into the irreducible 
representations of the SO (2) symmetry group results in 
summations over even indices only: 

C30 

S2{m= a2„(£)e^2"r (2) 

m— — oo 

Such a decomposition is deemed useful when each of the 
scalar functions aim (^) is itself a homogeneous function of 
its argument, characterized by an m dependent exponent: 

|a2m(A^)| ^ A^^'^|a2™WI , (3) 

where | • | stands for the norm of a complex number. For 
an isotropic fracture in an isotropic medium we expect 
a2m(^) — for all m ^ 0. In usual mode I experiments 
in which the crack propagates along the x direction and 
the tensile load is in the normal direction, there should be 
the same physics along lines with angles Q and ~d. This 
invariance under d — > —0 implies that the arguments of 
all aimi^) 7^ should be or tt. In reality this invari- 
ance might not hold when the fracture process breaks the 
symmetry dynamically; see bellow for an example. 

Example: aluminum alloy - Our first example was 
obtained 4] from a compact tension specimen made of 
7475 aluminum alloy first precracked in fatigue and then 
broken under tension in mode I. The raw fracture surface 
and the second order structure function computed from 
it are shown in Figs. [3 and El respectively. One sees the 
anisotropy of Sii^) with the naked eye. To quantitatively 
characterize this anisotropy, the structure function was 
decomposed as in Eq. (|2Jl. The log-log plots of ao{l), 
2|a2(^)| and 2|a4(£)| are exhibited in Fig. 2J By perform- 
ing linear fit of the relevant range in the log-log plots we 
find the following exponents 

C^^^= 1.32±0.07 ,C2^^= 1.45±0.08 ,cf ^= 2.1±0.1 . (4) 




FIG. 4: Log-log plots of the amplitudes ao{£), 2\a2{£)\ and 
2|a4(£)| vs. £ for the aluminum alloy. 
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The implication is that at smaller length-scales the 

(2) 

smaller exponent Q should be dominant and vice versa. 
Indeed, examining again the contour plot in Fig. |3| one 
observes that at small scales the contours tend to ellipses 
of smaller eccentricity, whereas at larger scale the con- 
tours are ellipses with increasing eccentricity. 

The crucial test of this approach is whether one can re- 
construct the structure function 5*2 (€,0) using the func- 
tional form of the irreducible representation and a mini- 
mal number of parameters. Indeed, at smaller values of £ 
the first two irreducible representations suffice. Writing 

S2{i, 0) « 8.30 f-^^ + 3.22 ^1-45 €05(26* + n) , (5) 

we compare in Fig|Sl the experimental data to Eq. (O 
for ^ = 5, 15 and 25/zm. The excellent fit is obvious. In 
fact, with four parameters (two amplitudes and two expo- 
nents) we can represent the structure function to within 
1% in norm as long as < 30/zm. For larger values of i 
the agreement decreases, and we need to employ the next 
irreducible representation. Adding 0.026 ^^'^ cos(40 -I- tt), 
we find the fit shown in Fig. for £ — 35fim. Beyond 
these values the power-laws fits lose their credibility for 
this experimental data set. 

Example: artificial rock- The second example was 
obtained from the dynamic fracture of artificial rocks pro- 
duced from carbonatic aggregates cemented by epoxy . 
The samples are plates of size 400 x 400 x 9 mm, and 
the fracture surface was measured using a scanning laser 
profilometer. The analysis of the experimental data fol- 
lows verbatim the first example. The plots of ao(^) and 
2|a2(^)| are shown in Fig. |3 Fitting the plots we find 

C^^^= 0.86±0.05 , C2^^= 1.93±0.05 , cf ^= 1.93±0.1 . (6) 

Two comments are in order. First, one should notice the 
non-universality of the scaling exponents as compared 
with the previous example This is to be expected 
when comparing a quasi-static and a rapid fracture ex- 
periments. Second, the present surface does not satisfy 
a, 6 ~9 symmetry. This results from the dynamic in- 
stability that leads to side branches that leave directed 
traces on the fracture surface, destroying the 9 —9 
symmetry 0. Due to the lack of symmetry the ampli- 
tudes of the coefficient am can take any phase, not con- 
strained to or TT as required by the 9 —9 symmetry. 
The lack of symmetry is clearly obvious in the reconstruc- 
tion of the structure function from the irreducible repre- 
sentations. We point out the advantage of the present 
approach in clearly revealing, through the lack of sym- 
metry, some aspects of the underlying physical process 
generating the fracture surface. In Fig. |S1 we compare 
the experimental values of S2{i, 9), for £ = 2 mm, to the 
expansion 

S2i£,9) « 0.025 0.0016 cos(26i -I- 2.09) 

+ 5.4 X 10"'' i^-^^cos{49 - 0.17) , (7) 
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FIG. 5: The experimental S2{£,d) for the aluminum alloy 
(circles) and the representation Eq. Q (line), for £ = 5, 15 
and 25/^m. 
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FIG. 6: The experimental S2{£,0) for the aluminum alloy 
(circles) and the SO (2) expansion up to the third even order 
irreducible representation (line), for £ — Sbfim. 
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FIG. 7: Log- log plots of the amplitudes ao{£) and 2|a2(^)| for 
the artificial rock. 
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FIG. 8: Comparison of the experimental S2{£,0) (circles) and 
the SO(2) expansion up to the third even irreducible repre- 
sentation (line) for the artificial rock with i = 2mm. 

The fit is satisfactory and the asymmetry in 9 is obvious. 

Taking the present two examples as representative, it 
appears that the S0(2) decomposition extracts pure scal- 
ing behavior in each sector, but that the scaUng expo- 
nents are not universal, at least in the two examples 
discussed here. Considering cuts in S2{£,9) along the 
9^0 and = tt directions, the present approach pre- 
dicts a mixture of scaling exponents rather than pure 
power laws. We verified that plotting the expansion 
for = and 6 = tt/2 in log- log plots we can reproduce 
the apparent scaling laws obtained by longitudinal and 
transverse cuts. Finally, one should point out that the 
SO (2) decomposition is not expected to yield satisfactory 
results when the material itself is strongly anisotropic. As 
an example we considered fracture surfaces in wood. This 



is clearly an anisotropic medium due to the fiber struc- 
ture, and indeed we found that along and across the fiber 
structure the scaling behavior appears credible, whereas 
the S0(2) decomposition failed altogether to reveal clean 
scaling properties in any sector. 

Summary: we propose that materials which can be 
fractured in an isotropic fashion, i.e materials having an 
isotropic structure, often have anisotropic fracture sur- 
faces only because of the breaking of isotropy by the ini- 
tial conditions. In such cases it appears useful to analyze 
the anisotropic contributions as "corrections to scaling" 
beyond the isotropic sector, which is always there, with 
a leading scaling exponent. The analysis reveals non- 
universality in the scaling exponents, a finding that calls 
for further future study and assessment, including the in- 
teresting question of the possible existence of universal- 
ity classes. On the practical side, we have demonstrated 
that the full information concerning the two dimensional 
structure function can be efficiently parameterized by a 
few amplitudes and scaling exponents. The reader should 
note that this Letter dealt only with second order struc- 
ture functions. In analogy to turbulence it may be possi- 
ble to decompose any higher order structure function into 
S0(2) irreducible representations 0. This may reveal 
additional interesting scaling properties such as the phe- 
nomenon of multiscaling (2| ■ Finally, we would like to em- 
phasize the great interest in the proposed isotropic frac- 
ture experiment and the measurement of the roughness 
exponent in such an experiment. If indeed this scaling 
exponent were identical to the exponent of the isotropic 
sector in a standard experiment, this would significantly 
strengthen the theoretical interest in the approach pro- 
posed in this Letter. 
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